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A possible ground state of NaV20s is proposed based on the Hartree approximation for both 
on-site and intersite Coulomb interactions. The results indicate that the intersite Coulomb 
interaction induces a zigzag type of charge disproportionation (i.e. charge ordering) along the 
ladders of V-ions resulting in the localized spins between neighboring ladders to form a spin 
gap. This new state, which is different from the spin-Peierls state so far believed, seems to be 
consistent with the existing experimental results. 
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Recently NaX^Os has been attracting interest. B 
Based on early structural study which postulated two dif- 
ferent V sites, i.e. V 4+ -ions and nonmagnetic V 5+ -ions,El ) 
this compound has been considered to be described by 
a simple spin 1/2 Heisenberg chain consisting of V 4+ re- 
sulting in the spin-Peierls (SP) state. 

However, there are many controversial issues as re- 
gards the nature of this material. The X-ray superlat- 
tice reflection below the transition temperature T c = 35 
K where the magnetic susceptibility steeply decreases, 
shows that the unit cell gets four times larger along the 
c-axis (the chains are along the 6-axis), which is not ex- 
pected from simple SP dimerized lattice distortion.!}* On 
the other hand, in an inelastic neutron scattering mea- 
surement below T c , two exitation branches along the a*- 
axis together with an intensity modulation were observed 
along the a* -axis, which can not be explained based on 
simple ID SP chainsif Furthermore, a recent crystal 
structural analysis showed that the V-ions are all equiv- 
alent at room temperature thus have valences 4.5+, BIS* 
while a 51 V-NMR experiment on single-crystalline sam- 
ple of NaV205 revealed that these V sites change its va- 
lences at T c from uniform V 4,5+ to two different sites of 
V 4+ and V 5+ .LP At the same time, A/T c , where A is the 
spin gap energy, is significantly large compared to other 
SP compounds, e.g. TTF-CuBDT, MEM-(TCNQ) 2 and 
CuGeOs.^ 

The causes of these puzzling experimental results is 
poorly understood, and even the origin of the observed 
insulating behavior in the whole temperature range is 
not clear yet. As regards the orbitals of <i-electrons at V 
sites, those with d xy symmetry are suggested to be rele- 
vant above and below T^Stl and then this d^-band is 
nominally 1/4-filled. There are theoretical studies based 
on the assumption of diaxge ordering into chains of V 4+ 
and V 5+ on one hand,BE£P while there is a proposal that 
each electron is located in a bonding V-O-V molecu- 
lar wavefunctionBliiP However, as iias-|been clarified in 
charge-transfer organic compounds ,t2l \}3 there exist var- 
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ious possibilities of charge distribution in 1/4-filled sys- 
tems triggered by the intersite Coulomb interaction when 
the geometry and consequently transfer integrals of the 
system are very anisotropic. 

In this letter, the possobility of a charge ordering in 
the ground state of NaX^Os has been pursued by use 
of the Hartree calculations including intersite Coulomb 
interaction. Furthermore, the origin of the spin gap be- 
havior is discussed based on these results. 

To investigate the effect of Coulomb interaction, the 
V2O5 layer, namely its a-b plane is considered, where the 
spatail arrangement of V-ions is schematically shown in 
Fig. 1. They form two- leg ladders in the 6-direction 
and these ladders align in the a-direction; a unit cell 
contains four vanadiums as in Fig. 1. The hopping pa- 
rameters between them, t a: tb : t xy: 
intersite Coulomb interactions, U and V% 
are taken into account; we define Vij = V a ,Vb,V a 



the on-site and the 
j, respectively, 



xy 



as 



the interaction along the bond with hopping t a ,tb,t xy , 
respectively, as in Fig. 1. Our Hamiltonian is a model 
of single orbital, d xyi for each vanadium site given as 
follows; 
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where < z, j > denotes the neighbor site pair, a is a spin 
index which takes j and j, and a\ a (a^) denote the 
number operator and the creation (annihilation) opera- 
tor for the electron of spin <j at the zth site, respectively, 
and rii — +riji. 

The Coulomb interactions U, Vij are treated within 
the Hartree approximation in a manner similar to that 
in refs. 12 and 13, 

n^nn -> (n iT ) riii + n *T ( n a) ~ ( n ^) ( n a) 
riifij — > (m) rij + rii (rij) - (n») (rij) , (2) 

and self-consistent solutions are searched for. Our calcu- 
lations are carried out at T = and the average electron 
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density is fixed at 1/4 per vanadium site. Several types of 
antiferromagnetic (AF) solutions are obtained and their 
energies are compared so that the true ground state is 
determined. The total energy £ is calculated as 

£ = rrr— V" eikaKF (Clka) 

K T ) (n zi ) - ViJ (m) (%•) , (3) 

i <hj> 

where N ce u is the total number of the cells, e^ a is the Zth 
eigenvalue at each k and rip (e) is the Fermi distribution 
function. 

The AF solutions we found to be stable are shown in 
Fig. 2. The 'dimer' type AF (Fig. 2(a) ) has uniform 
charge density, i.e. 0.5 electron per site, and, if the pairs 
of vanadium are taken as dimers, it can be considered as 
AF between these dimers along the ladders, which cor- 
responds to the state assumed in refs. 5 and 11 In the 
'chain' type AF (Fig. 2(b)), charge disproportionation 
(charge ordering) exists; the charges concentrate on one 
leg of the ladders, and AF ordering emerges along these 
chains. This type corresponds to the conventional ID 
V 4+ spin 1/2 chains. The 'zigzag' type AF is the one 
newly found, which also has the charge disproportion- 
ation but the sites with more electrons, i.e. V 4+ , are 
distributed in a zigzag way along the ladders as in Fig. 
2(c). In the actual solutions obtained by Hartree cal- 
culations, there exist, to be precise, finite but small spin 
moments on the sites which are described as open circles. 
However, as we will see below, within realistic range of 
parameters, these spin moments are small so that they 
can be neglected. 

For the explicit calculations the values of parameters 
are taken from ref. 11 as t a = 0.35 eV, tb = 0.15 eV, 
t X y = 0-3 eV and U = 4.0 eV, and Vij are varied main- 
taining the condition V = V a = Vb = V xy /V2. Here we 
assumed Vij oc 1/dij, where dij is the distance between 
the < z, j > pair, and approximated that the distances 
between intraladder vanadiums are all same and the dis- 
tance of interladder vanadiums is times smaller than 
that. The calculated energies of each state as a function 
of V for three different types of AF solutions together 
with that for the paramagnetic solution are shown in 
Fig. 3. This result shows that 'zigzag' type of AF has 
the lowest energy within a range of realistic values of V, 
namely V ~ 1 eV. The magnitudes of spin moment on 
sites 1 and 2 (see the inset of Fig. 4 where a unit mag- 
netic cell is shown), S z (i) (z=l,2), and the amount of 
charge disproportionation 5 in this case of 'zigzag' type 
AF are shown in Fig. 4. We note that the charge densi- 
ties are 0.5+5 on sites 1 and 4, 0.5—6 on sites 2 and 3, 
and the spin moments have the relation S^(l)= — S z (4) 
and S z (2)=— S z (3). When V is decreased, this 'zigzag' 
state gives away continuously to the 'dimer' state at 
V = V c , where S Z (1)=S Z (2)=-S Z (3)=-S Z (4). When 
V>0.5 eV, which is expected to be realistic, S ~0.5, i.e. 
the charge is almost concentrated in sites 1 and 4, and 
|S*(l)| = |Ss(4)| -1.0 (|S*(2)HS*(3)| ~ 0). Hence we 
can consider that the electrons are almost fully localized 
on sites 1 and 4 to form V 4+ while V 5+ in sites 2 and 3. 



If quantum fluctuation is taken into account in this state 
where localized 1/2 spins exist on the sites 1 and 4, it 
is easy to imagine that the pairs of nearest spins, which 
is along the bond with hopping t xy , form essentially a 
localized singlet. 

There exists a different proposal for the values of the 
hopping integrals with very small t xy , i.e. t a = 0.38 eV, 
t b = 0.17 eV and t xy < 0.1 eVB* We note, however, the 
result of our Hartree calculations even for t xy = eV 
does not show qualitative difference from those given in 
Figs. 3 and 4. We also note that the results are not 
largely affected by the condition V a = Vb = V xy /\/2 that 
we adopted; the results for V a = Vb = V xy show the same 
features. 

Now we compare our present results of the ground 
state with the 'zigzag' pattern of charge ordering, and the 
various experimental results noted above. The spin gap 
behavior in the magnetic susceptibility measurements 
can be understood as the spin singlet formation as dis- 
cussed above. TJbis is not a SP state and then the com- 
parison of A/T c u to other SP compounds is irrelevant. 
The change in- valences through T c observed in the NMR 
measurementEP may be due to the phase transition from 
the 'dimer' state above T c to the 'zigzag' state below T c . 
The 'zigzag' type charge ordering result in the period of 
two V sites in the 6-direction and of four ladders in the a- 
direction (see Fig. 2(c)), namely doubles the original unit 
cell size in the a-b plane Q which is consistent with the ob- 
served results below T c .u The exitation spectra from the 
'zigzag' singlet ground state is an unanswered question, 
but a deviation from ID character is expected. Thus it 
is natural that the inelastic neutron scattering measure- 
ment showed a dispersion which is not expected from ID 
SP chains. IS However, more experimetal studies, espe- 
cially the firm determination of crystal structure both 
below and above T c are needed for further discussions, 
e.g. regarding the unit cell size enlargement along the 
c-directiorJJ and the intensity modulation along the_al- 
axis in the inelastic neutron scattering measurement .urtS 

In summary, the ground state of NaV^Os is investi- 
gated within the Hartree approximation for both on-site 
and inter site Coulomb interactions. The effect of inter- 
site Coulomb interaction results in a zigzag pattern of 
V 4+ along the ladders of vanadiums resulting in the sin- 
glet pairs, which is different from the conventional pic- 
ture that this compound is a spin-Peierls system. This 
proposed ground state is consistent with the observed 
spin gap and also with the other experimental results so 
far reported. 
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Fig. 1. A schematic representation of the structure of V ions in 
the a-b plane of NaV^Os, where V ions are represented as black 
circles, together with the hopping integrals (after the notations 
of ref. 11), and the intersite Coulomb interactions Vij. Thick, 
thick-dotted, and thin-dotted lines representing the hopping pa- 
rameters ta^b and t X y, respectively. 



Fig. 2. A schematic view of different types of AF ordering for 
the structure given in Fig. 1. Circles and arrows represent V- 
ions and the spin moments, and +5 and —5 denote the amount 
of charge disproportionation (charge ordering) 0.5+5 and 0.5—5, 
namely more or less electrons than the average, respectively. 



Fig. 3. The energy per unit magnetic cell of each state shown in 
Fig. 2. The inset shows the region of small V. 



Fig. 4. The spin moment, S z , and the degree of charge dispro- 
portionation, 5, for the zigzag type AF as a function of V. The 
inset shows the magnetic unit cell along the ladder with the site 
indices, where the alignement of the spin moments are indicated 
by arrows and the amount of charge is 0.5+5 and 0.5-5, respec- 
tively. 



(a) 'dimer' 




(c) 'zigzag' 



